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The radial-orbit instability is a collective phenomenon that has heretofore only been observed in spherical systems. We 
find that this instability occurs also in triaxial systems, as we checked by performing extensive iV-body simulations whose 
initial conditions were obtained by sampling a self-consistent triaxial model of a cuspy galaxy composed of luminous 
and dark matter. iV-body simulations show a time evolution of the galaxy that is not due to the development of chaotic 
motions but, rather, to the collective instability induced by an excess of box-like orbits. The instability quickly transforms 
such models into a more prolate configuration, with 0.64 < b/a < 0.77 and 0.6 < c/a < 0.7 for the dark halo and 
0.64 < b/a < 0.77 and 0.59 < c/a < 0.67 for the luminous matter. Stable triaxial, cuspy galaxies with dark matter 
halos are obtained when the contribution of radially-biased orbits to the solution is reduced. These results constitute the 
first evidence of the radial-orbit instability in triaxial galaxy models. 
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1 Introduction 

The aim of the current work is to test, by extensive A-body 
simulations, the stationarity of the triaxial galaxy models re- 
alized by Capuzzo Dolcetta et al. (2007, hereafter CLMV07). 
These models were constructed by means of the standard 
orbital superposition method introduced by Schwarzschild 
(1979), and represent triaxial, cuspy elliptical galaxies em- 
bedded in triaxial dark halos in which the dark and the lu- 
minous components have the same axial ratios (a = 1, b = 
0.86, c = 0.7, i.e., maximal triaxiality). In CLMV07 these 
models were referred to as MODI and MODI -bis; the only 
difference between the two solutions was the maximum in- 
tegration time permitted during the orbital integration, which 
was longer in MODI -bis (5 Hubble times) than in MODI (2 
Hubble times), i.e., MODl-bis represents a more stationary 
solution. 



2 The Self-Consistent Triaxial Models 

The Schwarzschild method used to build the self-consistent 
solutions consisted in minimizing the discrepancy between 
the model cell "analytical" masses (obtained by integration 
of a given p[x, y, z)) and the masses given by a linear com- 
bination of orbits computed in the potential generated by p. 
In practice, the quantities to be, separately, minimized are : 
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where B 



k]j ;lm(dm) 



is the fraction of time that the fcth or- 
bit spends in the jth cell of the luminous-matter grid (dark- 
matter grid); Mj-i m is the mass which the model places in 
the jth cell of the luminous-matter grid and Mj- t d m is the 
same quantity for the dark-matter grid. Ck-im and Ck-dm 
represent the total mass, respectively, of luminous compo- 
nent and dark matter placed over the fcth orbit (1 < fc < 
n or b)- The basic constraints are Ck-im > and Ck-dm > 0, 
i.e., non-negative orbit weights. 

The mass model considered in CLMV07 for the lumi- 
nous component was a triaxial generalizations of Dehnen's 
spherical model (Dehnen 1993) with a weak cusp. The den- 
sity law for the luminous component was 
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and M the total luminous mass. 

For the dark component the adopted mass density was 
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and 
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dm^k^j;dm 



fc=l 



and pdmfl is the central dark matter density (Burkert 1995). 
Therefore the dark component has a flat, low-density core. 
In the present work all quantities are given in units corre- 
sponding to ai = Mi = t cross = 1 where t cross is the 
(2^ialf mass crossing time of the system. Assuming Mi = 
1O U M0, ai = 1 kpc and t cross = 26.5 Myr. 
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MODEL 


Solution 


L 


N lm 


N dm 


HL 


MODI 


23.71 


19684 


144886 


HLtis 


MODI - bis 


23.58 


19709 


146485 


LL 


MODI 


0.40 


19684 


144886 


LLbis 


MODI - bis 


0.34 


19709 


146485 



Table 1 Initial features of models used for the TV-body 



simulations. TV 



Irn(dm) 



is the number of particles used to 



represent the luminous(dark) component; L is the absolute 
value of the total angular momentum. 



3 N-body initial conditions 

In order to study the dynamical properties of MODI and 
MODI -bis we sampled these models as TV-body systems 
and let them evolve. The initial conditions for the TV-body 
integrations were set populating the fcth orbit with a number 
of particle proportional to Ck and randomly choosing posi- 
tions and velocities from the files containing the results of 
the orbital integration. 

In this work we define cc(z)-tube orbits as orbits having 
a non-vanishing x(z) component of the angular momentum. 
All other orbits (either box or chaotic) are classified as semi- 
radial orbits. Regarding the initial conditions of the TV-body 
systems, particular care should be given to the direction of 
rotation of the tube orbits. In fact the persistence of the sense 
of motion on tubes might be cause of internal streaming 
motions, giving a "rotating" model (Schwarzschild 1979; 
Merritt 1980). As a consequence of the symmetries of the 
potential, the orbits are invariant to a change in sign of the 
velocity and a decision should be made on the sense of mo- 
tion of the particles placed on tube orbits. To detect the ef- 
fects of any rotational motion, and to distinguish them from 
deformations due to dynamical instabilities, we set two JV- 
body simulations for each solutions(MODl and MODI -bis) 
in the two extreme cases of high and "zero" angular momen- 
tum. In Table Q] some parameters of four TV-body systems, 
sampling MODI and MODI -bis, are given. 

4 Results 

Simulations reveal that both MODI and MODI -bis evolve 
in shape, with no particular differences detected between 
the two cases during the evolution. The variation of the axial 
ratios is shown in Fig.Q] All the systems were found to have 
a final quasi-axisymmetric prolate shape. As example, Fig. 
|2] shows snapshots of the luminous component for model 
LLbi S . 

We evaluated the anisotropy parameter 2T r /T t where 
T r =< v 2 r /2 > and T t =< v 2 /2 >. For all TV-body sys- 
tems, roughlty the same values of the anisotropy parame- 
ters were found [2T r /T t ] dm « 2 and [2T r /T t ]i m w 1.4. 
These values are quite high, in particular for the dark com- 
ponent, suggesting that a bias in the semi-radial orbits may 
have induced an analog to the radial-orbit instability (ROI) 
of spherical models. 



5 Minimizing the contribution from 
semi-radial orbits 

To show that the dynamical instability seen in the TV-body 
simulations is due to the large radial "pressure" in the mod- 
els, new solutions containing a lower fraction of semi-radial 
orbits were computed by means of the orbital superposition 
technique and then evolved in time as TV-body systems. Our 
goal was to demonstrate the existence of a threshold value 
of the ratio between the number of tube orbits to that of 
semi-radial orbits above which stability is guaranteed. 

Following Poon and Merritt (2004), the relative orbital 
abundance in the self-consistent models was changed by 
adding a new term in equations ((TJ and (|2), which become 
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As above, the basic constraints were Ck- t i m > and Ck 
, i.e., non-negative orbit weights. 

Here, Wk-im{dm) is the "penalty" associated with the 
fcth orbit of the luminous (dark) component, having the ef- 
fect of filtering the orbital content in the solution: as Wk.j m rd m ) 
increases, the mass contribution Ci.-i m (d m ) of the fcth orbit 
in the model decreases. We chose Wk-i m = Wk;dm = for 
the tube orbits and W k -i m = W R -i m > and W k - dm = 
Wfi-dm > for the semi-radial orbits. The new solutions 
were found by means of the full set of orbits of the MOD1- 
bis catalog. Four different TV-body models were obtained 
by sampling new solutions, taking randomly the sense of 
rotation for the particles on tube orbits. Each of these mod- 
els was built using about 20000 luminous and 150000 dark 
matter particles. Table [2] gives the anisotropy of these TV- 
body systems: as expected, [2T r /T t ]i m and [2T r /T t ]dm are 
decreasing functions of Wr. 

6 The new simulations 

Fig. |3] shows the axis-length evolution of the new models. 
Model C is clearly dynamically unstable, evolving into a 



MODEL 


WR-lm 


WR-dm 


[2T r /T t ]lm 


[2T r /T t ]dm 


A 


50 


50 


0.512 


1.175 


B 


5 


5 


0.784 


1.335 


C 


5 x 10~ b 


5 x 10~ J 


1.220 


1.754 


D 


5 x 10~ b 


50 


1.230 


1.174 



Table 2 Anisotropy parameters of the new models. 
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Fig. 1 Upper panels: evolution of the axial ratios of /iX(left) and LL(right). Lower panels: evolution of the axial ratios of H Lua (left) and 
LL bis (right). R is the distance from the center where the axial ratios are evaluated. 
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Fig. 2 Time evolution of the luminous matter for model LL bia , pro- 
jected onto the x — y plane. 



is reduced. The correlation between the bar formation and 
the value of the radial velocity dispersion is a clear sign 
that the dynamical instability, seen above, can be identified 
with the ROI. Furthermore, the results of model D suggest 
that the instability disappears when the dark halo gets more 
isotropic, i.e., the ROI derives mainly by the initial high ra- 
dial pressure in the dark matter halo. 

Finally, we can draw the following conclusions from our 
analysis of the various different runs: 

i The dynamical instability that characterizes the time evo- 
lution of the galaxy models built in CLMV07 disappears 
when solutions contain a smaller number of semi-radial 
orbits; therefore it can be identified with the ROI; 

ii The ROI is due to the high concentration of radially "bi- 
ased" orbits in the dark matter halo. In particular nu- 
merical experiments suggest that stable configurations 
are obtained when [2T r /T t ]d m < 1-4. 



prolate configuration. By contrast, for models A, B and 
D, we did not observe the quick transition to instability 
that characterizes the evolution of models constructed in 
CLMV07. Therefore, we conclude that the instability can 
be suppressed when the radial contribution to the solution 



7 Discussion 

Our study has important implications for future construc- 
tion of self-consistent models of realistic galaxies. In par- 
ticular, it has been shown that the dynamical properties of 
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solutions cannot be directly deduced by the Schwarzschild Schwarzschild, M., 1979, ApJ, 232, 236 
method itself since the temporal development of models can 
be strongly affected by the distributions in velocity space. 
In this context, A^-body tests of the stability of the mod- 
els seems to be indispensable tool for drawing conclusions 
about general properties of models constructed via orbital 
superposition. 

Regarding the ROI, because a complete understanding 
of its mechanism is still lacking, it would be interesting 
to investigate more in depth its role in the case of triaxial 
potentials. Our simulations show dynamical features which 
cannot be detected in the case of spherical symmetry. The 
method exploited here to change the semi-radial orbit con- 
tribution to the models may be useful for constructing other 
radially-unstable triaxial systems. 
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